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Abstract 

The paper contains a simplified and improved version of the results obtained 
by the authors earlier. Wave propagation is discussed in a network of branched 
thin wave guides when the thickness vanishes and the wave guides shrink to a one 
dimensional graph. It is shown that asymptotically one can describe the propagating 
waves, the spectrum and the resolvent in terms of solutions of ordinary differential 
equations on the limiting graph. The vertices of the graph correspond to junctions 
of the wave guides. In order to determine the solutions of the ODE on the graph 
uniquely, one needs to know the gluing conditions (GC) on the vertices of the graph. 

Unlike other publications on this topic, we consider the situation when the spec- 
tral parameter is greater than the threshold, i.e., the propagation of waves is possible 
in cylindrical parts of the network. We show that the GC in this case can be ex- 
pressed in terms of the scattering matrices related to individual junctions. The 
results are extended to the values of the spectral parameter below the threshold 
and around it. 



1 Introduction 

Consider the stationary wave (Helmholtz) equation 

H^u = —e'^Au = Am, x E ^le, Bu = on dVt^^ (1) 

in a domain fig C i?*^, d > 2, with infinitely smooth boundary (for simplicity), which has 
the following structure: is a union of a finite number of cylinders Cj,s (which will be 
called channels) of lengths Ij, 1 < j < N, with diameters of cross- sect ions of order O {e) 
and domains Ji^^, ■ ■ • , JM,e (which will be called junctions) connecting the channels into 
a network. It is assumed that the junctions have diameters of the same order 0{e). The 
boundary condition has the form: B = 1 (the Dirichlet BC) or B = (the Neumann 
BC) or B = e-^ + a(x), where n is the exterior normal and the function a > is real 



*The authors were supported partially by the NSF grant DMS-0706928. 



1 



valued and does not depend on the longitudinal (parallel to the axis) coordinate on the 
boundary of the channels. One also can impose one type of BC on the lateral boundary 
of Qg and another BC on the free ends (which are not adjacent to a junction) of the 
channels. 

The axes of the channels form edges Tj, 1 < j < N, of the limiting {e 0) metric 
graph r. The junctions shrink to vertices of the graph F when e — > 0. We denote the set 
of vertices Vj by V. Let m channels have infinite length (m = for bounded Q^)- We 
start the numeration of Cj^^ with the infinite channels. So, Ij = oo for 1 < j < m. 




Figure 1: An example of a domain fl^ with four junctions, four unbounded channels and 
four bounded channels. 

The goal of the paper is the asymptotic analysis of the spectrum of if^, the resolvent 
{He — A)~^, and solutions of the corresponding non-stationary problems for the heat and 
wave equations as e — > 0. One can expect that -ffg is close (in some sense) to a one 
dimensional operator on the limiting graph F with appropriate gluing conditions (GC) at 
the vertices v & V. The ODE on F appears in a natural way from the following principle: 
the oscillating modes in the wave guides survive as e —>■ and the exponentially decaying 
and growing modes disappear. However, the justification of this fact is not always simple. 
In order to determine the solutions of ODE on F uniquely, one needs to know the GC 
on the vertices of F. The form of the GC in the general situation was discovered quite 
recently in our papers |2ii|-l23j. It turned out that they can be expressed in terms of 
scattering matrices for problems of the wave propagation through individual junctions of 
Qg. These GC hold in all the cases: in the bulk of the spectrum A > Aq, and near the 
threshold A ~ Aq, for bounded and unbounded fig. 

Equation ([T]) degenerates when e = 0. One could omit in ([T]). However, the problem 
under consideration would remain singular, since the domain shrinks to the graph F 
as £ — > 0. The presence of this coefficient in the equation is convenient, since it makes the 
spectrum less vulnerable to changes in e. As we shall see, in some important cases (spider 
domains Q^) the spectrum of the problem ([T]) does not depend on e, and the spectrum in 
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the same cases will be magnified by a factor of if e"^ in ([T]) is omitted. The operator 
= —e^A introduced in ([1]) will be considered as the operator in L'^{Q^). 
An important class of domains fls are the self-similar domains with only one junction 

and all the channels of infinite length. We shall call them spider domains. Thus, if Qs is 

a spider domain, then there exists a point x = x{e) and an e- independent domain Q such 

that 

= {{x + ex) : X G fi}, (2) 

i.e. Qe is the e-contraction of = 

For the sake of simplicity, we shall assume that is self-similar in a neighborhood of 
each junction. Namely, let Jj{v),e be the junction which corresponds to a vertex v & V oi 
the limiting graph V. Consider a junction ^ = Jj{v),e and all the channels adjacent to 
J^^s- If some of these channels have finite length, we extend them to infinity. We assume 
that, for each f G V, the resulting domain which consists of the junction J^^f, and the 
semi-infinite channels emanating from it, is a spider domain. We also assume here that all 
the channels Cj^e have the same cross-section u;^. This assumption is needed only to make 
the results more transparent (The general case is studied in [22]). From the self-similarity 
assumption it follows that is an e— homothety of a bounded domain uo C W^'^. 

Let Aq < Ai < A2... be eigenvalues of the negative Laplacian —Ad-i in uj with the 
BC Bqu = on duj where Bo coincides with the boundary operator B on the channels, 
see (lH), with e = 1 in the case of the third boundary condition. Let {fn{y)}, y & E 
R'^~^, be the set of corresponding orthonormal eigenfunctions. Then A„ are eigenvalues of 
— £:^Arf_i in uj^ and {e~'^^'^ipniy/£)} are the corresponding eigenfunctions. In the presence 
of infinite channels, the spectrum of the operator H,. consists of an absolutely continuous 
component which coincides with the semi-bounded interval [Ao,c>o) and a discrete set of 
eigenvalues. The eigenvalues can be located below Aq and can be embedded into the 
absolutely continuous spectrum. We will call the point A = Aq the threshold since it is the 
bottom of the absolutely continuous spectrum or (and) the first point of accumulation 
of the eigenvalues as e — 0. Let us consider two of the simplest examples: the Dirichlet 
problem in a half infinite cylinder and in a bounded cylinder of length /. In the first case, 
the spectrum of the negative Dirichlet Laplacian in fi^ is pure absolutely continuous and 
has multiplicity n + 1 on the interval [An, 00). In the second case the spectrum consists of 
the set of eigenvalues Xn,m = A„ + e^rn? /t^, n > 0, m > 1. 

It was shown in [2T]-[22] that the wave propagation governed by the operator if^, 
£ — s> 0, as well as the asymptotic behavior of the resolvent [H^ — X)~^ and of the eigen- 
values of ifg above Aq can be described in terms of the scattering solutions. While many 
particular cases of that problem with A = Aq + 0{e^) or A < Aq were considered (see 
[2]- [28]), the publications [21]- [23] were the first ones dealing with the case A > Aq, and 
the first ones where the significance of the scattering solutions for asymptotic analysis 
of was established. In particular, it was shown there that in both cases A > Aq and 
A ~ Aq, the GC of the operator on the limiting graph F will be expressed in terms of 
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the scattering matrices of the auxihary problems on the spider domains associated to 
individual junctions. A more profound analysis of the case A ~ Ag can be found in [23] . 

The main goal of this paper is to overview the results of [22]- [23] and simplify the 
proofs. We will mostly deal with the case of A G (Aq, Ai) where the results and proofs are 
more transparent. The number of scattering solutions is the smallest in this case and the 
scattering matrix is of the smallest size (compared to the case A > Ai). One of our main 
results is as follows. 

Theorem 1. // Aq < A < Ai, then the resolvent {He^^ — A)^^ can be approximated by 
{He^^ — (A — Ao))^\ where He^^ = —e^-^ is the operator of the second derivative defined 
on functions <; on the limiting graph T with the GC of the form 

ie[h + r,(A)]^^^(0) - v/A^[/, - T,(A)]^^(0) = 0, 

Here T^(A) is the scattering matrix of the axillary problem on the spider domain which 
corresponds to the junction J^^^, and (;'" is the vector which consists of restrictions of the 
function q (defined onT) onto edges adjacent to v. 

To be more exact, for any compactly supported f , the following relation is valid on 
channels outside of the support of f with exponential accuracy 

{H, - Xy'f ~ [{H^ - (A - Xo))-'fo]My/e), e^O, fo =< f, Vo{y/e) > . 

A more accurate statement of this theorem as well as some of its generalizations will 
be given in section 5. 

Note that the eigenvalues of the problem in are located not only below the threshold, 
but also above it. They depend on e and move very fast on the A-axis as e — > 0. Thus one 
can not expect to obtain an asymptotic approximation of the resolvent [H^ — A)^^ when 
A = A' > Aq is fixed and e —>■ 0. An asymptotic approximation of the resolvent (ifg — A)~^ 
as £ — > can be valid only if an exponentially small (in e), but depending on e, set on 
the A-axis is omitted. Another option is to fix A = A' > Aq and pass to the limit as e ^ 
without 6 taking values in some small set which depends on A'. 

While the condition A > Aq is natural for the wave propagation, the properties of the 
heat and diffusion processes depend on spectrum of Hs near A = Aq. As a by-product of the 
simpler approach to the problem introduced below, we will get a better result concerning 
the asymptotic behavior of the eigenvalues of if^ in bounded domains fl^ as e 0, 
A = Aq + 0(5^). It was shown in [22], [23] that the main terms of the eigenvalues of 
when A = Aq + 0{e'^), e — *• 0, coincide with the eigenvalues of the operator on the limiting 
graph with the GC defined by the scattering matrix at A = Aq. An explicit description of 
GC at A = Aq for arbitrary junctions (of order 0(e)) was also given there. Significantly 
later (see publications in arXiv), some of our results where repeated in [Hj. The new 
elements there are the location of the eigenvalues below the threshold and more accurate 
asymptotics of eigenvalues near the threshold. We will show here that the approach used 
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in [22] and [23] provides an approximation of the eigenvalues near the threshold with an 
exponential accuracy as well as the location of the eigenvalues below the threshold. 

The plan of the paper is the following. The elliptic problem in with a fixed e = 1 
will be studied in the next section. In particular, the scattering solutions are defined 
there. The asymptotic behavior of the resolvent {H^ — A)^^, of the spectrum and of the 
scattering solutions as £ — > 0, A > Aq, are obtained in section 3 for the simplest domains 
with one junction (spider domains). The one dimensional problem on the limiting graph 
will be studied in section 4. The case of arbitrary domains is considered in section 5. 
The last section is devoted studying the spectrum near the threshold. 

2 Scattering solutions and analytic properties of the 
resolvent when £ is fixed. 

We introduce Euclidean coordinates {t, y) in channels Cj^s chosen in such a way that 
the t-axis is parallel to the axis of the channel (so, t is not time, but space variable!), 
hyperplane Ry~^ is orthogonal to the axis, and Cj.e has the following form in the new 
coordinates: 

Cj^e = {it,ey) : < t < /j, ye u}. 

If a channel Cj^^ is bounded {Ij < oo), the direction of the t axis can be chosen arbitrarily 
(at least for now). If a channel is unbounded, then t = corresponds to its cross-section 
which is attached to the junction. 

Let us recall the definition of scattering solutions for the problem (fl]) in when 
A G (Ao,Ai). Consider the non- homogeneous problem 

{-6^ A -\)u = f, xe n,; Bu = on dn,. (3) 

Definition 2. Let f G L^^^IQ^) have a compact support, and A < Ai. ^4 solution u of 
^ is called outgoing if it has the following asymptotic behavior at infinity in each infinite 
channel Cj^^, 1 < J < 

u = ajc' ^^ipoiy/e) + 0(e"^), a > 0. (4) 

Remarks. 1. Here and everywhere below we assume that 

Im^A - Ao > 0. (5) 

Thus, outgoing solutions decay at infinity if A < Aq. 

2. Obviously, if (jlj) holds with some a > 0, then it holds with any a < VAi — A. 

Definition 3. Let A < Ai. A function = '^p \ 1 < p < m, is called a solution of the 
scattering problem in if 

{-6^ A - A)^ = 0, xe n^] = on dn„ (6) 



5 



and \1/ has the following asymptotic behavior in infinite channels Cj^ei 1 < J < ""^ : 

^l^^ = [6p,e-''^' + t,,e^'^']My/e) + 0{e-f), t ^ oo, a > 0. (7) 
Here 6pj is the Kronecker symbol, i.e. 6pj = 1 if p = j, Spj = if p ^ ]■ 

Remark. If Aq < A < Ai, then the term with the coefficient 5pj in ([7j) corresponds to the 
incident wave (coming through the channel Cp^e), the terms with coefficients tpj- describe 
the transmitted waves. The transmission coefficients tpj- = tpj{e, A) depend on e and A. 
The matrix 

T = [tp,] (8) 

is called the scattering matrix. Note that the scattering solution and scattering matrix 
are defined for all A < Ai. We assume that Im-y/A — Aq > when A < Aq, and the incident 
wave is growing (exponentially) at infinity in this case. 

The outgoing and scattering solutions are defined similarly when A G (A„, A„+i) (see 
|22j). In this case, any outgoing solution has n + 1 waves in each channel propagating 
to infinity with the frequencies \/X — Xg/e, < s < n. There are m(n + 1) scattering 
solutions: the incident wave may come through one of m infinite channels with one of 
(n + 1) possible frequencies. The scattering matrix has the size m{n + 1) x m{n + 1) in 
this case. 

Standard arguments based on the Green formula provide the following statement. 

Theorem 4. When Aq < A < Ai, the scattering matrix T is unitary and symmetric 
(tpj = tj^p). 

The operator iJ^ = — e^A, which corresponds to the eigenvalue problem ([1]), is non- 
negative, and therefore the resolvent 

Rx = {H, - A)-i : L\Q,) ^ L\Q,) (9) 

is analytic in the complex A-plane outside the positive semi- axis A > 0. If fi^ is bounded 
(all the channels are finite), then operator Rx is meromorphic in A with a discrete set 
A = of poles of first order at the eigenvalues A = /i^ ^ of operator H^. If fl^ has at 

least one infinite channel, then the spectrum of has more complicated structure (see 
Theorem [5] below). In this case, the spectrum has an absolutely continuous component 
[Aq, oo), and the resolvent Rx is meromorphic in A G C\[Ao, oo). We are going to consider 
the analytic extension of the operator Rx to the absolutely continuous spectrum. One can 
extend Rx analytically from above (ImA > 0) or below, if it is considered as an operator 
in the following spaces (with a smaller domain and a larger range): 

Rx:Ll^m^LUn,). (10) 

These extensions do not coincide on [Ao,C)o). To be specific, we always will consider 
extensions from the upper half plane (ImA > 0). We will call ffTOj) truncated resolvent of 
the operator H^, since it can be identified with the resolvent multiplied from the left 
and right by a cut-off function. 
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Theorem 5. Let have at least one infinite channel. Then 

(1) The spectrum of the operator Hg, consists of the absolutely continuous component 
[Ao, oo) and, possibly, a discrete set {iij^s} of non-negative eigenvalues A = /ij^e > with 
the only possible limiting point at infinity. The multiplicity of the a.c. spectrum changes 
at points A = A„, and is equal to m{n + 1) on the interval (An, A„+i). 

(2) The operator / fJ^) admits a meromorphic extension from the upper half plane 
ImX > onto [Aq, oo) with the branch points at X = of the second order and poles of 
first order at X = fij^,;. In particular, if A„ G {fJ'j,e}, then operator / flOj) has the form 



(3) If f & Llg^{Qe), A < Ai, and X is not a pole or the branch point (X = Xq) of the 
operator / fl^) . then the problem is uniquely solvable and the outgoing solution u 
can be found as the L^^^(fie) limit 

u = Rx+iof. (11) 

(4) There exist exactly m different scattering solutions for the values of X < Xi which 
are not a pole or the branch point of the operator fT^) . and the scattering solution is 
defined uniquely after the incident wave is chosen. 

(5) The scattering matrix T is analytic in A, when X < Ai, with a branch point of 
second order at X = Xq and without real poles. 

The matrix T is orthogonal if X < Xq. 

Remark. Let A„ ^ {/^i.s}- If the homogeneous problem with A = A„ has a nontrivial 
solution u such that 

u = ajifniv/e) + 0(e"'^*), x e Cj^s, t ^ oo, 1 < j < m, 7 > 0, (12) 

then Rx+io = f}'^\ +0(1), A — A„. If such a solution u does not exist, then operator 
f fTOj) is bounded in a neighborhood of A = A„. 

Proof of Theorem [5]. All the statements above concern the problem with a fixed value 
of e and can be proved using standard elliptic theory arguments. A detailed proof can be 
found in [23; a shorter version is given below. 

In order to prove the part of statement (1) concerning the absolutely continuous spec- 
trum of the operator H = —A, we split the domain fi^ into pieces by introducing cuts 
along the bases t = of all infinite channels. We denote the new (not connected) domain 
by fi'j, and denote the negative Dirichlet Laplacian in Q'^ by H'^, i. e. H'^ is obtained 
from Hi; by introducing additional Dirichlet boundary conditions on the cuts. Obviously, 
the operator H'^ has the absolutely continuous spectrum described in statement (1) of the 
theorem. Since the wave operators for the couple if^, H'^ exist and are complete (see [1]), 
the operator H,r has the same absolutely continuous spectrum as H'^. 

The remaining part of statement (1) and statements (2) and (3) can be proved by one 
of the well known equivalent approaches based on a reduction of the boundary problem 
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([3]) to a Fredholm equation which depends analytically on A. For this purpose one can use 
a parametrix (almost inverse operator): equation ([3]) is solved separately in channels and 
junctions, and then the parametrix can be constructed from those local inverse operators 
using a partition of unity (allowing one to glue the local inverse operators), see [22]. A 
similar approach is based on gluing together these local inverse operators using Dirichlet- 
to-Neumann maps on the cuts of the channels which were introduced in the previous 
paragraph. 

Statements (4) and (5) follow immediately from statement (3) and TheoremHl Indeed, 
one can look for the solution "^p^ of the scattering problem in the form 

^i'^=Xe-''^'My/e)+u (13) 

where x ^ C°°{Qs), X = outside Cp^e, X = 1 i^i C*p_e H {t > 1}. This reduces problem 
(E]), dZD to problem Q, (jlj) for u with / supported on Cp^ fl {0 < t < 1}. Statement 
(3) of the theorem, applied to the latter problem, justifies statement (4). Function u, 
defined in ( IT3I1 . satisfies the homogeneous equation ([3]) in infinite channels Cj^e, 3 7^ Vi 
and in Cp^^ fl {t > 1}, and it is meromorphic at the bottoms of these channels (at t = 
for j 7^ p, and t = 1 when j = p). Solving the problems in these channels by separation of 
variables, we obtain that the scattering matrix T is meromorphic in A, when A < Ai with 
a branch point of second order at A = Aq. It also follows from here that T is real valued 
when A < Aq. Analyticity of T and Theorem H] imply that T is orthogonal when A < Aq. 
From the orthogonality (A < Aq) and unitarity (Aq < A < Ai) of T it follows that T does 
not have poles. □ 



3 Spider domains, e ^ 0. 

Let us recall that is called a spider domain if it is self-similar (see ([2D) and consists of 
one junction and several semi-infinite channels. 

Theorem 6. Let be a spider domain and A < Ai. Then 

(1) the eigenvalues A = /ij.e = yUj of operator and the scattering matrix T do not 
depend on e, 

(2) the truncated resolvent / fJOj) has the following estimate: if f is supported on e- 
neighborhood of the junction, then 

\Rxf\ < C6-'s-'''/^\\f\\L2, A < Ai, S=dzst{\,{f^,}), (14) 

outside of 2e -neighborhood of the junction, 

(3) the scattering solutions have the following form on the channels of the domain: 

^(^) = [5p,,e-^'^* + tp,,e^'^Vo(l//^) + r^,,, a; G 1 < j < m, (15) 

where Ir^jl < Cd^^e^"^^ when £ > 0, | > 1 and < A < Ai. Here a < VAi — A is 
arbitrary, C = C{a). 
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Remark. Formula (fT5|) looks similar to the definition ([7]). In fact, the remainder in 
(j7]) decays only when t —* oo, and ([7j) does not allow us to single out the main term of 
asymptotics of scattering solutions as e — > 0. 

Proof. All the statements above follow immediately from self-similarity of the domain 
fls- Namely, we make the substitution 

x^^^ (16) 

(see (E])) and reduce problem ([3]) in to the problem in Q which corresponds to e = 1. 
These two problems have the same eigenvalues and scattering matrices. This justifies the 
first statement. Let Vx, g be functions R\f,f after the change of variables (fT6|) . From 
statement (2) of Theorem [5] it follows that 

WvxWlHk) < C6-'\\g\\L^. = CrV^/2||/||i2, 



where K consists of the parts of the channels of Q where 1 < t < 3. Then the standard a 
priori estimates for the solutions of the equation Au — Xu = imply the same estimate 
for \vx\ on the cross sections t = 2 : 

<CrV^/2||/||r2, t = 2. 



\vx\ S e 

The latter implies the same estimate for If^l when t > 2 by solving the equation Au — Xu = 
in the corresponding regions of the channels of Q with the boundary condition at t = 2. 
This justifies the second statement of Theorem [61 The last statement can be proved 
absolutely similarly. We reduce the scattering problem in to the scattering problem in 
Q and use representation ( |T3l) with e = 1. This implies ((71) with e = 1 and the remainder 
term rpj such that \rpj\ < C6^^e~°'^ for t > 1. It remains only to make the substitution 
inverse to (fT6l) . □ 

In spite of its simplicity. Theorem [61 allows us to obtain two very important results: 
small e asymptotics of the spectrum of and the resolvent {H^ — A)~^ for arbitrary 
networks of thin wave guides Q^- For this purpose, we need to rewrite (1TB]) in a slightly 
different (less explicit) form. 

We denote by <jpj the linear combination of exponents in the square brackets in f[T5l) . 
This is a function on the edge Tj of the graph. Let be the vector-column with compo- 
nents (^pj, 1 < j < m. Obviously, ^p satisfies the equation 

(^'^ + A - Ao)^ = 0. (17) 

We will use notation "^^^ for both the scattering solution and the column-vector whose 
components "^^j are restrictions of the scattering solution "^p^ on the channels Cj,£, 1 < 
j < m. Then (1151) can be rewritten in the vector form as 



= '^pMy/e) + r(^) = [e,e-^^* + tpe^'^']My / e) + (18) 
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where x G Ui<j<mC*j,e, is the vector with components r^j, all components Cpj of the 
vector Cp are zeroes except Cp^p which is equal to one, and tp is the p-th column of the 
scattering matrix T. Let us construct the mxm matrix with columns \E'p'^^ and the matrix 
<^ with columns (^p, I < p < m. As it is easy to see, ^(0) = (/ + T), ^-'(0) = i^^^^f^{-I + T), 

and therefore, 

te{I + r)^'(O) - - Ao(/ - T)^(O) = 0. (19) 

Of course, this equality also holds for individual columns <jp of matrix 

It is essential for extending the results to arbitrary networks of wave guides that 
the gluing condition f|T9|) together with some condition at infinity is equivalent to the 
explicit form of qp given by f|T8|) . Namely, let ^ satisfy f|T7|) . Then 

^ = OpC * - * + /3pe' = 

with some constant vectors Op, /5p. We will say that q = tpp is a solution of the scattering 
problem on the graph T with the incident wave coming through the edge Tp if ipp satisfies 
equation f|T71) . GC f|T9l) . and Op = Cp, i.e.. 



= epe-'^* + /3pe^^* (20) 

Thus, we specify the incident wave and impose the GC defined by the scattering problem 
in Q^, but we do not specify the scattering coefficients of the outgoing wave. The next 
theorem shows that the scattering problem on the graph will have the same scattering 
coefficients as the problem on fi^. 

Theorem 7. Formulas l[TS\) . ( fl^) and "^^^ = ippipoly/e) + Vp'' are equivalent. 

Proof. It was already shown that qp defined in f fTSj) satisfies f|T9|) . Conversely, if we write 
(3p in fl20|) as tp + hp and put fl20|) into f|T9l) . we will get that hp = 0, i.e. ipp coincides with 
<^p defined in ([H]). □ 



4 One-dimensional problem on the graph. 

The spectrum of the operator Hg, an the asymptotic behavior of the resolvent will be 
expressed in terms of the solutions of a problem on the limiting graph F which is studied 
in this section. 

Let fie be an arbitrary (bounded or unbounded) domain described in the introduction, 
and let F be the corresponding limiting graph. Points of F will be denoted by 7 with t 
being a parameter on each edge Tj of the graph. We are going to introduce a special 
spectral problem 

•= -^'^^ = - ^0)^ (21) 
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on smooth functions <j = <j(7) on F which satisfy the following GC at vertices. We split 
the set V of vertices v of the graph into two subsets V = V1UV2, where the vertices from 
the set Vi have degree 1 and correspond to the free ends of the channels, and the vertices 
from the set V2 have degree at least two and correspond to the junctions Jy^e- We keep 
the same BC a.t v E Vi as at the free end of the corresponding channel of fie (see (II])): 

B^ = atv e Vi. (22) 

The GC at each vertex v E V2 will be defined in terms of an auxiliary scattering 
problem for a spider domain fly^- This domain is formed by the individual junction J^^g 
which corresponds to the vertex v, and all channels with an end at this junction, where 
the channels are extended to infinity if they have a finite length. Let T = T„(A) be the 
scattering matrix for the problem ([T]) in the spider domain Q'^^^ and let ly be the unit 
matrix of the same size as the size of T. We choose the parametrization on F in such a 
way that t = at f for all edges adjacent to this particular vertex. Let d = d{v) > 2 
be the order (the number of adjacent edges) of the vertex v E V2. For any function ^ on 
F, we form a column- vector (t) with d{v) components which is formed by the 

restrictions of ^ on the edges of F adjacent to v. We will need this vector only for small 
values of t > 0. The components of the vector q^""^ are taken in the same order as the 
order of channels of Q'^ ^. The GC at the vertex f G V2 has the form 

d 



ie[Iy + Ty{X)]-&\t) - - Ao[/. - T,(A)]^('^)(t) = 0, t = 0, veV2, (23) 

if A 7^ Aq. Condition fl23|) can degenerate if A = Aq, and it requires some regularization in 
this case. 

Solutions of fl2Tl) have the following form 



= aje' ^ * + bje ' ^ , 7 e ^j- 

If ImA > and G L'^iT), then bj = for infinite edges (see ([5])). Thus, if satisfies 
equation (l2Ti) in a neighborhood of infinity, then 



<^ = aje'^^^\ 7 G Tj, I < j < m, t » 1. (24) 

We will assume that condition ( 12^ holds also when A is real, i.e., we consider only those 
solutions of (I2T]) with real A = A' > Aq which can be obtained as the limit of solutions 
with complex X = X' + ie when e ^ 0. 

We will call function g = gx{'y,^;e), 7,^ G F, the Green function of the problem 
(21\)-(24\) if it satisfies the equation (with respect to variable 7) : 

-e':^g-iX-Xo)g = 5^{^), (25) 

and conditions fl22l) - fl24l) . Here ^ is a point of F which is not a vertex, and ^^(7) is the 
delta function supported on 7 = ^. 
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Lemma 8. Let A < Ai, A 7^ Aq. Operator /i^ = —e^-^ is symmetric on the space of 
smooth, compactly supported functions on T which satisfy conditions [2^) and [23\) . 
Proof. One needs only to show that 



^^,rit),,rit)j - = o, t = o, ve v,, (26) 

for any two vector functions <j = c; = c:^^^ which satisfy GC (1231) (similar relation at 
V G Vi obviously holds). Let A G (Aq, Ai). Then matrix Tt,(A) is unitary (Theorem Hj). If 
matrix + Ty is non-degenerate, we rewrite ([23]) in the form ji&\t) = Ac^^^'^t), t = 0, 
where the matrix 

A = ^^7 + T,(A)]-i[4 - T,(A)] 



is real. The latter immediately implies ( 1261) . Similar arguments can be used if 1^ — Ty 
is non-degenerate. If both matrices are degenerate (i.e., Ty has both eigenvalues, ±1), 
we consider a unitary matrix U such that UTyU* is a diagonal unitary matrix. Since 
{17^1,1/^2) = (<^i,<^2) for any two vectors ^1,^2, one can easily reduce the proof of fl26|) to 
the case when Ty is a diagonal unitary matrix. Then fl23|) implies the following relations 
for coordinates qj{t) of the vector (^^""^t) : gj(0) = 0^9(0) or ^^(0) = 6j?"j(0), where 
constants real. The first case occurs if the corresponding diagonal element of Ty 

differs from —1, the second relation is valid if this element is —1. These relations for qj{t) 
imply (1261) . Similar arguments can be used to prove (!26|) when A < Aq, since matrix Ty is 
orthogonal in this case (see Theorem [5]). □ 

Theorem 9. For any e > there is a discrete set A(e) on the interval [— Aq, Ai) such 
that the Green function gxi'y,^; e) exists for all A < Ai, A ^ ^i^): O'^'^d has the form 

- D{\e) ' ^^^^ 

where function h is continuous on the set 7,^ G F, A < Ai, e > and uniformly hounded 
on each hounded suhset, and 

D{\,e) = J:Z=icM)e''^'-- (28) 

Here Sm are constants, functions Cm(A) are analytic in X < Xi with a hranch point of 
second order at X = Aq, and D if X < Xq. 

Proof. We fix the parametrization on each edge Tj of the graph. Then, obviously, 

g, = a,e-''^' + b^e''^\ ^er^, if ^ ^ T,, (29) 

,. = a,e-^* + 6,e^^* + ^=4=sin[^^^(t-r)_], if ^ G F,. (30) 

V A — Aq ^ 
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Here r is the coordinate of the point ^, (t — r)_ = mm(t — r, 0), and the last term in (130!) 
is a particular solution of (l25l) on with a bounded support. There are 2N unknown 
constants in the formulas above where N is the total number of edges of the graph. 
Conditions ( |22i) - (l2^ provide 2N linear equations for these constants. As it is easy to see, 

the coefficients for unknowns in all the equations have the form a(A)e* ^ where a(A) 
is analytic in A < Ai with a branch point of second order at A = Aq, and s = or ±lj (Ij 
are the lengths of the finite channels). The exponential factors in the coefficients appear 
when the formulas fl29l) . fl30|) are substituted into GC at the end point of the edge Vj where 
t = Ij. We apply Cramer's rule to solve this system of 2N equations. This immediately 
provides all the statements of the theorem with D{\,e) being the determinant of the 
system. One only needs to show that -D 7^ for A < Aq. Note that the latter fact implies 
the discreteness of the set K{e) = {A : D{\,e) = 0}. 

Obviously, D = if and only if the homogeneous problem (12T]) - (12^ has a non-trivial 
solution. Let A < Aq. Then solutions ? of the problem decay at infinity, and 

= ^[-£V'-(A-Ao)?]?rf7 

= -e'K U + lie'Wr - (A - Ao)^^]ci7- (31) 

It was shown in the proof of Lemma [8] that it is enough to consider only diagonal 
matrices T when the terms under the sign S„ above with v E V2 are evaluated. Since T 
is orthogonal when A < Aq, the diagonal elements of T are equal to ±1. Then ( !23l) means 
that each component of the vector 

^(^) or its derivative is zero at the vertex. Hence, the 
terms in the sum above with f G V2 are equal to zero. They are zeroes also for those 
f G Vi where the boundary condition in fl22|) is the Dirichlet or Neumann condition. If 
f G Vi and B = e-^ + a, a > 0, these terms are non-positive. Hence, relation (13T]) implies 
that <^ = when A < Aq. □ 
Theorem [9] does not contain a statement concerning the structure of the discrete 
set A.{e). This set becomes more and more dense when e ^ 0. In general, every point 
A' g(Ao, Ai) belongs to A{e) for some sequence of e = £j(A') 0. However, it is not an 
absolutely arbitrary discrete set, but the set of zeroes of a specific analytic function fl28l) . 
and this fact provides the following restriction on the set A{e). 

Lemma 10. For each bounded interval [a, Xi], each o" > and some M, there are ce~^ 
intervals Ij of length a such that 

\D{X,e)\ > ca^' when e > 0, A G [a, Ai]\ U /,-, c = c{a). 

This Lemma is a particular case of Lemma 15 from [22] (the set Tq is empty in the 
case which is considered in this paper). 

In order to construct the resolvent of the problem in Q^, we need to represent the 
Green function gx of the problem fl25|) . fl22|) - fl24|) on the graph T through the solutions of 
the scattering problems on the spider subgraphs of T. 
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We will call a function ip = ippil) solution of the scattering problem on the graph T if it 
satisfies the equation (12T]) , conditions (l22l) , (123!) and has the following form at unbounded 
edges of the graph: 



M7) = Sp,je' - ' + apje' - *, 7 e T^-, l<p,j<m, (32) 

where Spj is the Kronecker symbol. This scattering solution corresponds to the wave 
coming through the edge Tp. These scattering solutions on the graphs were introduced in 
the previous section in the case when the graph corresponds to a spider domain. In fact, 
only this simple case will be needed below. 

Lemma 11. If the graph T corresponds to a spider domain fi^, then the scattering solution 
ippi'y) exists and is defined uniquely for all X < Xi, X ^ Xq. Any function <; on T which 
satisfies equation l[21\) and GC condition (E^ is a linear combination of the scattering 
solutions i'pi'j). 

Remark. For arbitrary graphs, one may have nontrivial solutions of the homogeneous 
problem (!2Tl) -( l24l) supported on the set of bounded edges of the graph. This occurs when 
A G A(£:). The set A{e) is empty for spider graphs. Proof. If we take apj = tpj, where 
tpj are the scattering coefficients in the spider domain Q^, then function fl32|) will satisfy 
fl23|l (see the derivation of f|T9l) ). Hence, the scattering solutions V'pIt) exist for all A < Ai, 
A 7^ Ao, since the scattering coefficients are defined for those A by Theorem [H If we put 
function ( 132|) with apj = tpj + hpj into GC ( l23l) . we immediately get that hpj = (see 
the proof of Theorem [7]). Thus, scattering solutions are defined uniquely. The space of 
solutions of equation (!2T]) is 2m-dimensional. The m x 2m dimensional matrix (/^ + T^(A), 
ly + Ty{X)) formed from coefficients in GC (l23l) has rank m. Hence, the solution space of 
the problem (1211) . (123!) is m-dimensional. Obviously, functions ipp are linearly independent 
on r. Thus any solutions of ( |2T]) . ( |23l) is a linear combination of functions ipp. □ 
Let Fjg be the edge of F which contains the point ^ (see (l25l)). We cut the graph F 
into simple graphs F(w) with one vertex v by cutting all the bounded edges at some points 
C,j G Tj. We will choose = ^. Let us denote by F'(t>) the spider graph which is obtained 
by extending all the edges of F(f ) to infinity. Let ipp^vil) be the scattering solutions on 
the graph F'(f ). 

Lemma 12. There exist functions 

a = ap^^{X,6,^), A < Ai, e > 0, ^ G Tj^, 
which are continuous, bounded on each bounded set, and such that 

ctp,?j(A, 5, ^) , , , , 

9^ = ^p—F^\ — ^'^pAV^ 7 e T{v). 
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Proof. It follows from the previous lemma that gx can be represented as a linear combi- 
nation of the scattering solutions: 

In order to find the coefficients Cp^y, we note that gx is equal to a combination of two 
exponents on the edge Tp C T{v) with the coefficient of the incident wave equal to Cp „ : 

= Cp,^e*^^* + bp^^e-'^^', jeTpC T{v). 
Now Cp^y can be found by comparing the formula above and (1271) at two points of Tp. □ 

5 Small £ asymptotics for the problem in Q^. 

As everywhere above, the domain Q^, considered below, can be bounded or unbounded. 
Denote by A° the union of eigenvalues of the operator ([3]) in all the spider domains fl'y ^ 
associated to ^2;.. These spider domains consist of individual junctions and all the channels 
adjacent to this junction. The channels are extended to infinity if they have a finite length. 
The set A° does not depend on e due to Theorem O Let us recall that A{e) is the set of 
eigenvalues of the one dimensional problem fl2Tl) - fl2^ on the limiting graph (see Theorem 

ED. 

The eigenvalues of the operator = —e'^A of the problem ([1]) which are located on 
the interval {—oo, Ai) are exponentially close to the set A''uA(e). In the process of proving 
this statement, we will get the asymptotic approximation of the resolvent {H^ — A)~^ as 

6 —>■ 0. Namely, the following theorem will be proved. 

Let A' < Ai and let A'^ be e^^-neighborhood of the set A° U A{e). Assume that the 
right-hand side of ([3]) has a compact support which is separated from junctions, i.e. there 
exist T,d > such that support of / belongs to UA^ where Aj is the part of the channel 
Cj £ defined by the inequalities r < t < Ij — t if Ij < oo, or t < t < d if the channel is 
infinite. 

Theorem 13. (1) There exists u > such that the eigenvalues fij^g of the operator 
which belong to the interval (— oo. A') with an arbitrary X' < Ai are located in e~J^- 
neighborhood of the set A° U A(£:). Here a = \i — A'. 

(2) Let the support of f belong to UAj and let u = Rxf be the solution of problem 
Here Rx is the truncated resolvent Then for any rj > 0, there exist u > and 
p = p{r]) > such that u = Rxf has the following asymptotic behavior in all the channels 
outside rj— neighborhood of the support of f 

u = Rxf=igxfo)Ml)+0{e^), A e (-oo, A')\A^ e^O. (33) 

Here 

/o = /o(7) = (/,^-'/Vo(f)), 7er, 
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and 'g\ is the integral operator on the graph T whose kernel is the Green function gx 
constructed in Theorem\^ 



9\fo = j^9x{l,^;£)fo{OdC 



Remark. Below, we also will get the asymptotics of u = R\f on the support of /, as 
well as a more precise estimate of the remainder in ( 133|) . 
Proof. Let 

fi = fiix) = f-e-'/'foM^), a;efi„ 

i.e. /o = /o(7) is the first Fourier coefficient of the expansion of / with respect to the 
basis {e~'^^'^(fj{-)}, and /i is the sum of all the terms of the expansion without the first 
one. We are going to show that u = R\f has the following form on the channels of fi^ : 

« = /?A/ = (?A/o)(/5o(|)+X^A/i + 0(e^), AG (-oo,A')\A^ £ ^ 0, (34) 

where j/, p > 0, x ^ C^iVLs) is a cut off function such that x = on all the junctions, 
X = 1 outside of e-neighborhood of junctions, and function -R°/i is defined by solving 
the following simple problem in the infinite cylinder. Let fij be the restriction of /i onto 
the channel Cj.e. We extend the channel Cj^e to infinity (in both directions) and extend 
fij by zero. Let Uj be the outgoing solution of the equation 

— e^Aw — \u = fij 

in the extended channel. Then R^fi is defined as R^fi = Uj in the channel C^.e Obviously, 
xRxfi can be considered as a function on Q^. 

The justification of (JMl) and the proof of the Theorem [13] are based on an appropriate 
choice of the parametrix ("almost inverse operator"): 

which is defined as follows 

Pxf = iGJo)M-) + ixRlfi) - ^.XvRUxvlie'A + X){xR'xfi) - fi]]- (35) 



Here L^^ is a subspace of L^ifle) which consists of functions supported on UAj. Now we 
are going to define and study, successively, each of the terms in the formula above. In 
particular, we need to show that 

-{e'A + X)Pxf = f + Qxf, Qx:Ll,-^Ll,, ||gA||<C7e^. (36) 

Operator Gx is an integral operator with the kernel Gx{x, z; e), x^z E Vt^, which is 
defined as follows. We split fi^ onto domains VL^^^ by cutting all the finite channels Gj^e 
using the cross-sections t = tj. Let z G Aj^. Then we choose tj^ to be equal to the 
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coordinate t = t{z) of the point z. Other cross-sections are chosen with the only condition 
that T < tj < Ij — r, i.e., the cross-section t = tj is strictly inside of A^. Let ^ be the 
spider domain which we get by extending all the finite channels of Qv,e to infinity. Let 
"^^l be the scattering solutions of the problem in the spider domain fi^^. The small e 
asymptotics of these solutions is given by Theorems [6] and [71 We introduce the following 
functions '^^Jj by modifying the remainder terms in these asymptotics: 

^ii = ^pMy/^) + Xv4'\ (37) 

where Xv ^ C°°{^e), Xv = ^ on r-neighborhood of the junction, Xv = outside of flv,e- 
Then we define G\ by the formula 

G,ix, z; e) = E.^^i^T^^S- ^ ^ ^''-^ (38) 
L>[A, e) 

where ap,„, D are defined in Lemmas fT2| [91 ^ is the point on the graph F which corresponds 
z G Ajp, i.e., the point on the edge Tj^ where t = tj^y Since function "^p}; satisfies the 
equation (e^A + X)u = on Qv,e, from Theorems [61 and [71 it follows that 

-(e^A + A)^^) =0(rVF), 5^0, -oo<X<X' , xen,,„ 
where a = Xi — X'. We choose < ^. Then 5 > for A G {—oo, X')\ A'", and 

-{e^A + X)¥;l = 0{e-'-^), e^O, X e {-oo, X')\ A" , x e n,,s- 

Since coefficients bounded. Lemma [TOl with a = e implies that 

- (£2A + A)G'a = 0(e-^), e->0, A G (-oo, A')\ A^ x G (39) 

Relations f[39l) are valid on each domain fit,,^. Now we are going to combine them and 
evaluate (e^A + X)Gx for all x G fi^. From ^ST\\ . f[551) and Lemma [T^ it follows that the 
function 

is infinitely smooth in the channels of Vt^. Here 7 is the point on F which corresponds to 
X G l^g. Then from (l39l) it follows that 

- (e^A + A)Ga = 5e(7)</'o(-) + 0(e-^), e ^ 0, A G (-00, A')\ A^ x G (40) 

As it is easy to see, the remainder in fHUl) is zero in the region where '^Xv 7^ 0, i.e., support 
of the remainder belongs to UAj. 

Now let us study the second and third terms in the left hand side of ( l35ll . Obviously, 

-{e'/\ + X){xRlh) = xh + h = h + K h= -2e'Vx-VRlh-e\Ax)Rlfi. (41) 
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Here we used the fact that x = 1 on the support of /i. Since /i is orthogonal to fo{-), 
function R^fi and all its derivatives decay exponentially in each channel Cj^s as - oo 
where r is the distance from Aj. Hence, 

h = 0ie-'^) = 0ie-'f), e^O, AG (-00, A'). (42) 

The remainder terms will be parts of the operator Qx, and we need the kernel of this 
operator to be supported on UAj. Unfortunately, h is supported on e- neighborhoods of 
the junctions. The last term in ( l35l) is designed to correct this. Since h is supported on 
the region where Vx 7^ 0, function h can be represented as the sum h = T^yh^, where 
hy = Xvh has estimate (H2|) and is supported on the ^-neighborhood of the junction Jy ^ 
which corresponds to the vertex v. Consider h = T,yXvRx vlXvh] which is defined as follows. 
We apply the resolvent i?° ^ of the problem in the spider domain Q'^ ^ to h^, multiply the 
result by Xv and extend the product by zero on fls\^v,e- 

From (H2I) and Theorem [6] it follows that 

|^A,AI ^ CS-^e-"^ < Ce"^, 5^0, AG (-00, A')\ A^ (43) 

if we choose z/ < |, so that S > . From standard a priory estimates for the solutions 
of homogeneous equation (e^A + X)u = it follows that estimate (HM is valid also for all 
derivatives of R\f, since this function satisfies the homogeneous equation outside of 2e- 
neighborhood of the junction. Then (H3l) holds for the derivatives of R^^h^. This allows 
us to obtain, similarly to fj^Tj) . that 

-{e^A + X)h = h + hi, /ii = 0(e-^), e 0, A G (-cx), A')\ A^ (44) 

where hi is supported on the closure of the set Vxv 7^ 0. This set belongs to UAj. Finally, 
from dinD, (BH), dBD it follows that 

-(e'A + \)Pxf = f + g, 9 = 0(6-^), e ^ 0, A G (-00, A')\ A^ (45) 

and g is supported on UA^. One can easily check that g depends linearly on /. Besides, 
one can specify the dependance on the norm of / in estimates of all the remainders above. 
This will lead to fl36l) instead of fH5l) . In fact, fl36|l is valid when Qx is considered as an 
operator in or as an operator in the space of continuous functions on UAj. 

We are going to construct now the solution u of problem Q with f E L'^^. We look 
for u in the form u = Pxg with unknown g E L"^^. Obviously, u satisfies the boundary 
conditions and appropriate conditions at infinity. Equation ([3]) in leads to g + Qxg = /• 
Since the norm of operator Qx is exponentially small, function g exists, is unique and 
9 = f + q, ||g||<Ce-f 11/11, i.e., 

u = Pxif + q), \\q\\Ll^<Ce-^\f\\Li^^, e ^ 0, A G (-00, A')\ 

This justifies ( l34l) and ( l33l) . The first statement of Theorem [T3] follows from here. Namely, 
assume that an eigenvalue /i = fij^s of the operator belongs to (— 00, A')\ A*^. Then 
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the truncated resolvent R\ (see ([9])) has a pole there (see Theorem [5]) . The residue of 
this pole is the orthogonal projection on the eigenspace of H^. The pole of R\f may 
disappear only if / is orthogonal to the eigenspace which corresponds to the eigenvalue 
X = fi. Non-trivial solutions of the equation (A + X)u = in fi^ can not be equal to zero 
in a subdomain of Q^- Thus, there is a function f E L'^^ which is not orthogonal to the 
eigenspace, and Rxf must have a pole at A = /i. This contradicts flM|) and fl33|) . □ 
The following statement can be easily proved using Theorem [13] and reduction f|T3|) of 
the scattering problem to problem ([3]), (jl]). 

Theorem 14. For any interval [a, X'),there exist p, z/ > such that scattering solutions 
\E'p^e(x) of the problem in Qe have the following asymptotic behavior on the channels of 
^le as e — > 

where ippi'y) = ipp^\'y) are the scattering solutions of the problem on the graph T and 

\4'\x)\<Ce^, Xe[a,X')\A\ 

Here 7 = 7(0;) is the point on T which is defined by the cross-section of the channel 
through the point x, and d{^) is the distance between 7 and the closest vertex of the graph. 



6 Eigenvalues near the threshold 

In some cases, in particular when the parabolic problem is studied, the lower part of the 
spectrum of the operator Hg, is of a particular importance. Theorem [T3] provides a full 
description of the location of the eigenvalues. One may have a finite number of eigenvalues 
below Aq. They are determined by the junctions and situated in an exponentially small 
neighborhood of A" (set of eigenvalues of the corresponding spider domains). The eigen- 
values between Aq and Ai are situated in an exponentially small neighborhood of K{e) 
which is the set of eigenvalues of the one dimensional problem (!2T|) - (p^ on the limiting 
graph (see Theorem [H]). 

We will assume that VL^ has at least one bounded channel (for example, VL^ is bounded). 
The opposite case is studied in Theorem [61 We also assume that A = Xo + Oi^e^). Then the 
eigenvalues of the problem f[2T]) - f[2^ will depend on the form of the GC f[23|) at A = Aq. 
Let put A = Aq + in f[2T]) - f[2^ . Then this problem takes the form 

- = on r, (46) 

5^ = aXve Vi, (47) 
^[/, + T,(Ao + p£')]^?('^)(t)-p[/.-T,(Ao + /i£')]^(''Hi) = 0, t = 0, v^V2, (48) 
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<^ = aje'>'\ 7 e Tj, l<j<m, t » 1. (49) 

The last condition is not needed if fi^ is bounded (m = 0). 

Since matrix T^(Ao) is orthogonal and its eigenvalues are ±1, the GC fj48l) with e = 
has the form 

P^W(0)=0, p^^^(-)(o) =0, veV2, (50) 

where P,P-^ are projections onto eigenspaces of matrix T^,(Ao) with the eigenvalues ^1, 
respectively. Let k be the dimension of the operator P, and d — k he the dimension of the 
operator P-*-, where d = d{v) is the size of the vector <^(''). Then (1501) imposes k Dirichlet 
conditions and d — k Neumann conditions on the components of vector c;^^^ written in 
the eigenbasis of the matrix T^^Xq). Note that the standard Kirchhoff conditions (^ is 
continues on F, a linear combination of derivatives is zero at each vertex) has the same 
nature, and k = d — 1 in this case. 

Problem fj46l) - fj49l) with e = has a discrete spectrum {/Xj}, j > 1, and the same 
problem with e > is its analytic perturbation. Thus, the following statement is valid. 

Theorem 15. If eigenvalues {/Uj} are simple, then eigenvalues {fij{6)} of problem [J^l- 
( [T^P are analytic in e : 

(^) = Yl ^J>^''' ^i'O = (51) 

n>0 

Remarks 1. This statement implies that eigenvalues A G A(e) in 0(£^)-neighborhood of 
Ao have the form 

A = A,(e) = Ao + £'^/i,>£". (52) 

n>0 

2. The assumption on simplicity of fij often can be omitted. For example fl5T|) . fl52l) remain 
valid without this assumption if k = d (the limiting problem is the Dirichlet problem). 
In the latter case one may have multiple eigenvalues (for example, when the graph has 
edges of multiple lengths), but the problem with e = is split into separate problems on 
individual edges. 

Theorem [T5] makes it important to specify the value of k in the condition (!50|) . This 
value depends essentially on the type of the boundary conditions at dfl^ and on whether 
A = Ao is a pole of the truncated resolvent ffTOl) or not. 

Definition 16. A ground state of the operator in a domain o,t X = Aq is the 

function ipo = 'ipoix), which is bounded, strictly positive inside Q^, satisfies the equation 
(— A — Ao) ^/'o = in (^iT-d the boundary condition on dQ^^ o.iT'd has the following 
asymptotic behavior at infinity 

ipo (x) = v^o (I) [pj + o (1)], X G Cj, \x\ +00, (53) 

where pj > and (fo is the ground state of the operator in the cross-sections of the 
channels. 
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Obviously, if the Neumann boundary condition is imposed on dQs, then Aq = 0, and 
the ground state at A = exists and equal to a constant. It was shown in [22], [23] that 
the ground state at A = Aq does not exist for generic domains fls in the case of other 
boundary conditions on dfli;. In particular, it does not exist if there are eigenvalues of 

below Ao or if the truncated resolvent does not have a pole at A = Aq. The following 
result was proved in [22], |23j . 

Theorem 17. (1) The ground state at X = implies k = d — 1. Thus the eigenvalues 
fj^jie), e 0, converge to the eigenvalues of the Kirchhojf problem in the case of Neumann 
condition on dQ,, (Q^ is arbitrary) and in the case of other boundary conditions on dQ,, 
for special, non-generic fl^. 

(2) If Dirichlet or Robin condition is imposed on dQ^ ond the truncated resolvent does 
not have a pole at X = Xo (this is a generic condition on Q^), then k = d and fij{e), 
e ^ 0, converge to the eigenvalues of the Dirichlet problem. 

Other possible (non-generic) GC at A = Aq are given by ( l50l) . 
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